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ABSTRACT 

Context. After the work of Gladman et al. (1998), it is now assessed that many irregular satellites are orbiting around 
Uranus. 

Aims. Despite many studies have been performed in past years, very few is know for the light-curves of these objects 
and inconsistencies are present between colours derived by different authors. This situation motivated our effort to 
improve both the knowledge of colours and light curves. 

Methods. We present and discuss time series observations of Sycorax, Prospero, Stephano, Setebos and Trinculo, five 
faint irregular satellites of Uranus, carried out at VLT, ESO Paranal (Chile) in the nights between 29 and 30 July, 2005 
and 25 and 30 November, 2005. 

Results. We derive light curves for Sycorax and Prospero and colours for all of these these bodies. 
Conclusions. For Sycorax we obtain colours B-V =0.839 ± 0.014, V-R = 0.531 ± 0.005 and a light curve which is 
suggestive of a periodical variation with period ft; 3.6 hours and amplitude ~ 0.067 ± 0.004 mag. The periods and 
colours we derive for Sycorax are in agreement with our previous determination in 1999 using NTT. We derive also a 
light-curve for Prospero which suggests an amplitude of about 0.2 mag and a periodicity of about 4 hours. However, the 
sparseness of our data, prevents a more precise characterization of the light-curves, and we can not determine wether 
they are one-peaked or two-peaked. Hence, these periods and amplitudes have to be considered preliminary estimates. 
As for Setebos, Stephano and Trinculo the present data do not allow to derive any unambiguous periodicity, despite 
Setebos displays a significant variability with amplitude about as large as that of Prospero. Colours for Prospero, 
Setebos, Stephano and Trinculo are in marginal agreement with the literature. 

Key words, planets and satellites: individual (Sycorax, Prospero, Stephano, Setebos, Sycorax, Trinculo), planets and 
satellites: general, methods: observational, methods: data analysis, methods: statistical methods: numerical, 

1. Introduction 

In recent years many irregular satellites has been discovered around Uranus (jGladman et al. I . Il998l : iKavelaars et al. I . 
2004; Gl adman et al. l . [2000tlSheppard. Jewitt. Klevna I . [2005Q . Irregular satellites are those planetary satellites on highly 
elliptic and/or highly inclined (even retrograde) orbits with large semi-major axis. These objects cannot have formed 
by circumplanetary accretion as the regular satellites but they are likely products of captures from heliocentric orbits , 
probably in association w i th th e planet formation itself (|Greenberg I . Il976t iMorrison fc BurnsI 119761 : iMorrison et al. I . 
Il977t I Jewitt fc Sheppardl . 20051). It is pos s ible for an object circling about the Sun to be temporarily trapped by a 
planet ( Heppenheimer I . Il975t iGreenbergl . Il976t IMorrison fc BurnsI . Il976t to cite only some). But to turn a tem- 



porary capture into a permanent one requires a source of dissipation o f orbital energy and that partic les could re- 
main inside the Hill sphere long enough for the capture to be effective dPollack. Burns fc Tauber" . 1979)) . Oth e rwise , 



the trapped object will escape within at most few hundred o f orbits (|Bvl fc Ovenden , Heppenheimer , 1975; 



Heppenh eimer fc Porco 1 . 119771 : iPollack. Burns fc Tauberl I1979D . During the planet formation epoch several mechanisms 
may have operated, some of which h ave the potential to be activ e even after this ear ly epoch. They fall mainly int o 
few categories: collisional interactions ([Colombo fc Franklin I . fl97ll ) . pull-down capture f Heppenh eimer fc Porco I . ll977j ). 



* Based on observations with the ESO Very Large Telescope + FORS2 at the Paranal Observatory, Chile, under program 
075.C-0023. 
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gas drag (jPollack. Burns fc Tauberl . Il979f ). four bodies interactions in the reference fra me of the the Sun-Planet system 
either between the captured body and a large regular satellite of the planet (|Tsui I . l2000f ). or between the two components 
of a binary object leading t o an exchange reaction wh ere one of the components of the binary is captured and the other 
is ejected from the system ([Agnor fe Hamilton! l2006f ). 

Collisional capture, the so called break-up process leads to the formation of dynamical groupings. The resulting 
fragments of the progenitor body after a break-up will form a population of irregular satellites expected to have similar 
composition, i.e., similar colours, and irregular surfaces. Large temporal variations in the brightness of irregular satellites 
are expected from rotating bodies of highly elongated shapes and/or irregular surfaces consistent with a collision fragment 
origin. 

Gas drag is expected t o occur in the environment of the proto planetary nebula (|Bvl fe Ovenden 1 . 119751 : iHoredt 1 . 119761 
Hcp penheimer fc Porco I . Il977t iPollack. Burns fc Tauber I . Il979f) and it may origin dynamical families of fragments. In 
this case fragments would be produced by the hydrodynam ical breaking of the intruding b ody in smaller chunks in 
case they exceed the tensile strength of the entering body (|Pollack. Burns fc Tauberl . |1979| ). Gravitational attraction 
of fragments prevents them from escaping, in general the hydrodynamical pull being not larger than self-gravity, but a 
small impact of a ~ 1 Km size object, like ly common in the nebula enviro nment, would be sufficient to disperse them 
without introduce a further fragmentation (|Pollack. Burns fc Tauberl . Il979f h A specific prediction of this scenario is the 
production of fragments with a more regular/round surfa ce than in the break-up process, leading to light-curve with 
low amplitude variations ( jPollack. Burns fc Tauberl . Il979t ) . 

On the contrary, if pull-down capture, four bodies interactions or exchange reactions are the dominant causes of 
formation of the irregular satellites, each object would be the result of an independent capture event. In this case, no 
obvious correlation between dynamical properties, colours and light-curves would be expected. 

To cast light on these scenarios, colours and light curves are very important, since they would allow one to discriminate 
between collisional or non-collisional origin for irregular satellites. Theories of irregular satellite capture have lacked 
many constraints. However, the rapidly-growing number of known irregular satellites is now providing new insights on 
the processes of planet formation and satellite capture. 

A possible origin of that the large obliquity of Uranus is a giant impact event between the planet and an Earth-size d 
planetesimal, occurred at the end of the epoch of accretion (|Slatterv. Benz fc Cameron! Il992t iParisi fc Brunini"! [1997) . 
The dynamical and physical properties of the Uranian irregular satellites may shed light on the ir capture mechanism 
and may also witn ess the mechanism leading to the peculiar tilt of the planet's rotation axis (jBrunini et al. I . I2002t 
IParisi et al. I . [20071 ). For example, significant fluctuations have been observed in the Caliban light-curve for which data 
are consistent with a light-curve with amplitude Acaiiban = 0.12±0.01 mag and a most probable period of » 2.7 hours as 
in th e Sycorax light c urve, As VC orax = 0-03 2 ± 0.008 mag with either a period of « 4.1 hours or w 3.7 hours ([Maris et al. I . 
1200 1|) . In this regard IRomon et al. I (|2001[ ) report discrepancies in the spectrum they possibly attributed to rotational 
effec ts. All of this se e ms to support the idea of a coll i siona l scenario. However, the existence of a dynamical group- 
in g (iKavelaars et al. I . l2004t ISheppard, Jewitt. Klevnal . [2005) is still debated on the light of the colour determination 
of iGrav et al. I ( 2004h . Regrettably, th ere seems to be a l ack of consistency between B-V and V-R colours of differ- 
ent authors for Sycorax and C aliban ([Maris et al. I . 120011 ; iRettig. Walsh fc Consolmagno I . 120011 : IRomon et al. I . [2001; 
IGrav. Holman fc Fraser I . l200l . This may be ascribed to systematic differences in the photometry and accompanying 
calibration or, at least for Caliban, to rotational effects. 

In an attempt to improve on the situation, in this paper we present and discuss new observations of five irregulars 
of Uranus, Sycorax, Prospero, Setebos, Stephano and Trinculo, obtained with the ESO Very Large Telescope on Cerro 
Paranal, Chile. 

The paper is organised as follow: Sect. [2] describes observations and data reduction, light-curves are discussed in Sect. 3, 
while Sect. 4 present the satellites colours. The conclusions are reported in Sect. 



2. Observations and Data Reduction 

We observed the i r regul ar satellites Sycorax, Prospero, Stephano, Trinculo and Setebos with the FORS2 camera 
(jAppenzeller et aT~l . ll998l ) at the focus of VLT Antu telescope in Paranal, Chile, in the two consecutive nights of July 28 
and 29, 2005. 

We used the standard FORS2 B, V, R, I filters Q, which are very close to the Bessel system. In particular, the 
effective wavelength, A c ff, and FWHM, A A, for the filters reported by ESO are A e ff,B = 0.429 m/i, AAb = 0.0880 m/z 
for B; A c ff,v = 0.554 m/i, AAy = 0.1115 m/i for V; A c ff3 = 0.655 m/i, AAr, = 0.1650 m/i for R; and X e s,i — 0.768 mfi, 
AAi = 0.1380 m/i for I. Stephano and Trinculo were observed in the same frames so that we observe five objects with 
just four sequences. For this reason Stephano Rl, R2, . . ., VI, V2, . . ., frames corresponds to Trinculo Rl, R2, . . ., VI, 
V2, . . ., frames. Each object has been observed in consecutive sequences of frames. After the end of the sequence for 
a given object the telescope switched to the sequence of another object. Ideally, colours would have to be calculated 
by combining magnitudes from frames in the same sequence, in order to limit the rotational effects. For each sequence, 
pointing of the telescope and orientation of the camera have been kept fixed. During the acquisition of each frame the 
telescope has been tracked at the same rate of the target, while the telescope has been resetted at the default pointing at 
the beginning of each frame in the sequence. However, given the slow proper motion and the short exposures, the effect of 
differential tracking on background stars has been negligible, background stars does not appear elongated. Interruptions 



1 See 



http://www.eso.org/instruments/fors/inst/Filters/ for further details. 
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due to a ToO, mid- night calibrations and some minor problem prevent us from keeping the same sequences in the two 
nights. Both nights have been photometric, with average seeing « 1.1 arcsec. FORS2 is equipped with a mosaic of two 
2k x 4k MI T CCDs (pixel si ze of 15 x 15 micron) with a pixel scale, with the default 2-by-2 binning, of 0" .2b /pixel. The 
satellite and lLandolt I (|1992f ) standard stars were centred in CCD #1. Pre-processing of images, which includes bias and 
flat field corrections, were done using standard IRAF0 routines. Aperture photometry was then extracted using the IRAF 
tool QPHOT, both for the standard stars and the satellite, using a handful of bright field stars to estimate the aperture 
correction. The resulting corrections were small, going from 0.06 to 0.25 in all filters. A series of R exposures have been 
taken with the aim to construct a light curve and search for some periodicity. A few B, V and / exposures have been 
taken as well to constrain the satellites' colours. The calibration was derived from a grand total of 3 standard stars 
per night in the PG0231+051, SA92, PG2331+055, MARK A, SAlll, PG1528+062 and PG1 133+099 |LandoltJ (fl992h 
fields. The two nights showed identical photometric conditions, and therefore a single photometric solution was derived 
for the whole observing run 

M = 77i inst + a m ~ m ■ Airmass; (1) 

where mi ns t = b, v, r or i are the instrumental magnitudes, M = B, V, R or J are the reduced magnitudes, Airmass is 
the airmass and a m and /3 m are the calibration coefficients. We obtain a m = 2.332, 2.864, 3.112, 2.546 and /3 m = 0.269, 
0.177, 0.147. 0.150 respectively for the B, V, R and I bands. No colour correction have been applied due to the very 
small colour term. A few additional observations of Prospero were acquired on the nights of November 22 and 25, 2005 
in compensation to the ToO. We reduced the data in the same way as in the July run, but obtained an independent 
calibration, being a m = 2.318, 2.879, 3.007, 2.529 for the B, V, R and I bands, respectively, which is very similar to the 
July one. The list of measures for the five satellites is in Tab. Q] The table reports the reduced magnitudes, errors and 
exposure times. The shortest exposures have been acquired to improve frame centering nevertheless we report magnitudes 
from these frames too. The time column refers to the starting time of each exposure. No corrections for light travel times 
have been applied to these data. 

3. Light Curves 

Fig.s[Tl [21 El Hand \5\ presents R light curves respectively for Sycorax, Prospero, Setebos, Stephano and Trinculo for the 
July 2005 nights. Each plot is splitted in two subpanels, the left for the July 29 th and the right for July 30 th . Squares 
in gray represents the measurements of a common field star with similar magnitude. To avoid confusion, error bars for 
the field star are not reported and the averaged magnitude is shifted. For the same reason the few data obtained in 
November 2005 for Prospero listed in the Tab[T]are not plotted in Fig. O 

We analyse magnitude fluctuations in the light-curves trying to assess first of all whether the detected variations 
may be ascribed to random errors, Sect. 13.11 instabilities in the zero point of the calibration, Sect. 13.21 or to the effect 
of unresolved background objects Sect. 13.31 In case the variability has been judged significant attempt to estimate the 
amplitude and the period, Sect. 13.41 We use both analytical methods and a Monte Carlo (MC) - Bootstrap technique. 

3.1. Testing against random fluctuations 

The results of a \ 2 test performed on V, R (and eventually B, I) measures are reported in the first two columns of 
Tab. [2] (similar results are obtained by a bootstrap on the data). In this test the hypothesis Ho to be disproved is that 
the data are compatible with a constant signal (different from filter to filter) with random errors as the sole cause of 
brightness fluctuations. The table shows that this hypothesis may be discarded for Prospero and Sycorax with a very 
high level of confidence. As for Stephano and Setebos the level is lower but still significant, whereas for Trinculo the 
hypothesis can not be discarded at all. Before considering the case for a periodic variation, the case for a systematic trend 
in the brightness is considered, since the irregular sampling in time prevents the application of robust de-trendization 
techniques. As evident from the table, even this case can be excluded by the present data at a level of confidence similar 
to the constant case. 

3.2. Field stars analysis 

In order to asses the level of calibration accuracy in an independent manner, several field stars having magnitudes 
encompassing those of t he irregular satellites and common to each frame in both nights have been measured in the same 
way as the satellites (see lCarraro et aL~l . l2006l for an example of the adopted technique). As shown in Fig. s [HO [21 HI and [5] 
field stars are characterised by less wide fluctuations than satellites. A more quantitative test is obtained taking a set of 
variability indicators by which measuring the variability of field stars and satellites and comparing them. This is done in 
Fig.[H]where two variability indicators, the peak-to-peak variation (top) and the RMS (bottom), for satellites (red spots) 
and field stars (light blue asterisks) are plotted against the R magnitude of the objects. The first important thing to 
note is the good consistency between the two indicators. While it is evident that the variability of Sycorax, Prospero and 
Setebos is above the level of variability of field stars, the same is not true for Stephano and Trinculo. In addition, test 
for the correlation of, as an example Sycorax or Prospero and the related field stars shows that they are not significantly 



2 IRAF is distributed by NOAO, which are operated by AURA under cooperative agreement with the NSF. 
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correlated. As an example the correlation coefficient between Sycorax and three field stars is Cs yC orax.s = —0.36, —0.08 
and —0.53. The probability that this level of correlation can be reached by chance even if their time series are not 
correlated are respectively 27%, 78%, 15%, while for Prospero Cp rospcr0iS = —0.10, 0.22 and —0.12 with probabilities 
respectively of 64%, 36% and 59%. In addition, even assuming a correlation between field stars fluctuations and satellite 
fluctuations, it would explain only a small fraction of the satellites variability. Indeed, variances of field stars accounts 
for 0.02% — 6% of the Sycorax variance, 6% — 9% of the Prospero variance, 12% — 18% of the Setebos variance, 2% — 14% 
of Stephano or Trinculo variances. In conclusion, field stars variability, connected to calibration instabilities, can not 
account for most of the variability of at least Sycorax, Prospero and Setebos, while as evident from Fig. [5] Stephano and 
Trinculo would have to be considered more cautiously. 

3.3. Unresolved background objects 

Unresolved objects in background may affect photometry. A quicklook to the framesEI shows that in general, the satellites 
passes far enough from background objects to allow a proper separation of their figures. It can be considered also the case 
in which a satellite crosses the figure of an undetected faint object in the background producing a fake time variability 
of the satellite light-curve. It is quite simple to compute the upper limit for the magnitude of a background object, 
-Rbck, able to produce the observed variations of magnitude for these satellites. The result is i?bck < 23.6 mag, 24.7 mag, 
25 mag, 25.5 mag and 25.5 mag, respectively for Sycorax, Prospero, Setebos, Stephano and Trinculo. Those magnitudes 
are within the detection limit for Sycorax and Prospero, near the detection limit for Setebos and marginally outside the 
detection limit for Stephano and Trinculo. In conclusion at least for Sycorax and Prospero, we can be confident that 
background is not important. 

3.4. Looking for amplitudes and periodicities 

Given the sparseness of our data it is not easy to asses safely the shapes, the amplitudes and the periods of our light- 
curves, despite at least for Sycorax, Prospero and probably Setebos a significant variability is present. However, we think 
important to attempt a recover of such informations, at least as a step toward planning of more accurate observations. 

Constraints on the amplitude for the part of light curves sampled by our data may be derived from the analysis of the 
peak-to-peak variation. After excluding data with exposure times below 100 sec we evaluate peak-to-peak variations for 
R, Ai? P 2p. Denoting with the amplitude of the light-curve we may assume A R ~ Ai? P 2 P /2. To cope with the random 
noise, Ai? P 2 P have been evaluated by MonteCarlo, simulating the process of Ai? P 2 P evaluation assuming that the random 
errors of the selected data are normally distributed. We obtain ^R.sycorax ?=i 0.07±0.01 mag, ^R.prosporo ii 0.27±0.04 mag 
and Ap^setebos ^ 0.31 ± 0.05 mag. Where the ^ symbol is used because the Ar w Ai? P 2 P /2 relation is strictly valid only 
if the true minima and maxima of light-curves are sampled, a condition which we are not safe to have fullfilled. Another 
order of magnitude estimate of the amplitudes is based on the analysis of their RMS, std(i? t ). It easy to realize that for 
any periodical light-curve of amplitude ^4r then std(i? t ) w /Dump Ar. Where /Dump > is a factor which depends both on 
the sampling function and the shape of signal. Assuming a sinusoidal signal sufficiently well sampled, /Dump ~ 1-39— 1.46, 
giving for Sycorax A R , Sy corax ~ 0.06 mag, while for Prospero and Setebos A RiProspcro « A R , Sc tcbos ~ 0.2 mag. 

We then consider the case for a periodical variation in our light-curves by attempting first to search for the presence 
of periodicities in the hypothesys of a sinusoidal time dependence, and in case of a positive answer, assessing the most 
likely sinusoidal amplitude. To cope with the limited amount of data increasing both the sensitivity to weak variations 
and the discrimination power against different periods, we fit the same sinusoidal dependence on V, R (and eventually 
B and I) assuming that colours are not affected by any significant rotation effect. In short the model to be fitted is 



where / — R, V (and eventually B and I) indicates the filter, Mf the measured magnitudes for that filter, M/\o the 
averaged magnitude for the filter /, A is the amplitudes (f> the phase and r an arbitrary origin in time. It has to be noted 
that phase-angle effects are not considered here. The reason is that no safe dependence of the magnitude on the phase- 
angle has been established so far for these bodies. On the other hand, the variation of the phase angles over two consecutive 
nights is just about 2.4 arcmin. Consequently, we expect the phase-angle effect to be quite negligible. As widely discussed 
in literature the problem of searching for periodicities by fitting a model with a sinusoidal time dependence is equivalent 
to the analysis of the periodogram for the given data set (|Lomb l . [l976tlScargle l . ll98"2tlCumming. Marcv fc Butler I . fl999; 
ICummingl , 12004 just to cite some). The Lomb and Scargle (hereinafter LS) periodogram is the most used version. In 
this view, a better formulation of the problem is obtained expressing the model in the following form: 



Mf = Acos(2ir/P(t - r) + <p) + M/, , 



(2) 




(3) 



3 A sequence of the observed frames for Prospero has been already published in iParisi et al. I (|2007T ). 
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where lu — ■jf, A c and A s co-sinusoidal and sinusoidal amplitudes related to A and (j> by the simple relations A = 

- y/ A 2 + A 2 and phas e <j> = arctan(A c /A s ). The time origin r can be arbitrarily fixed, but the canonical choice is 
(lLomb I . [19761 : iScargle I . [1981) 

l^f=B,V,B.,I 5^- 

tan(2wr) = = . (4) 

which will be our definition of r. Then, the free parameters involved in the minimisation are P, Mf^ for / =V, R (B,I), 
A c and A s or equivalently, A and </>. However, being interested to amplitude and not to the phase we marginalise our 
statistics over 4>. 

A set of best possible combinations of Mf,o for / =V, R (B,I), A c and ^4 S for each given P in a suitable range 
[•Pmin, Pmax] is obtained by minimising 

x 2 {p)= J2 cos (?(^ - T )) + ±&<t£ - T)) + MM - w ^) 2 (5) 

f=B,V,R,Ij=l f>i 

where TUfj are the magnitudes for filter / measured at the times tfj with associated errors er/j with j = 1, 2, . . ., Nf 
the index of Nf measures obtained for the filter /. The minimisation is carried out analy tically, with the r defined in 



Eq. |p. This way the method becomes a generalisation of the floating average periodogram (jCumming. Marcv fc Butlerl 
1999: Cummingt I2004D 



and reduces to it in the case in which magnitudes come from a single filter, while in the case of 
homoscedastic data and null zero points we return to the classical LS periodogram. 

The sensitivity to noise as a function of P is not constant and varies with time. Fig.[7]represents the result for a Montecarlo 
designed to asses the sensitivity to noise of the periodogram in wide range of periods for Sycorax (upper frame) and 
Prospero (lower frame). The Montecarlo code generates simulated time series assuming the same time sampling of data, 
the same errors, normal distribution of errors and time independent expectations, and then computes the corresponding 
periodogram. Dots in the figure represents the realization of such periodograms, the full-line represents the Fourier 
transform of the time window. The sparseness of that causes a strong aliasing with diurnal periodicities. The plot is 
dominated by the prominent 24 hrs diurnal peak, followed by the 12 hrs, 8 hrs, 6 hrs, 4 hrs, 2 hrs peaks of decreasing 
amplitude. It is evident that above a period of 20 hrs the sensitivity of the periodogram to random errors is rapidly 
increasing. Then, our data set is best suited to detect periods below 20 hrs, or better, due to the presence of the 12 hrs 
peak, periods below 10 hrs. 

Significant periodicities for P < 10 hrs are likely present in the data if at least for one P in the range, the squared 
amplitude A 2 (P) = A C 2 {P) + A 2 (P) obtained by minimising Eq. ©, exceeds a critical value A 2 rit , which is fixed by 
determining the false alarm probability for the given [P m i n ,Pmax], Pfa (^4 2 > ^crit) S- This interval has been sampled 
uniformly with a step size AP — 0.025 hours (the results do not depend much on the choice of the step size) and the 
Pfa as a function of ^4 2 rit has been assessed. In determining the Pfa we exploit the fact that we want to calculate this 
probability for values of A 2 rit for which Pfa is small. In this case, 

P FA (A 2 > A 2 cl , t ) » Pfa,o exp(-A c 2 rit /A c 2 rit!0 ), (6) 

which is good for Pfa < 0.3, and with the parameters Pfa,o and A 2 vit Q determined from Montecarlo simulations. The 
last two columns of Tab. [2] report the results of this generalised version of the LS periodogram. Again, for Sycorax and 
Prospero a quite significant periodical signal is detected. For Setebos the detection is marginal, while for Stephano and 
Trinculo no detection can be claimed at all. There are many reasons for the lack of detections of periodicity despite 
random noise cannot account for the variability. Among them a lack in sensitivity, the fact that the period is outside the 
optimal search window, and that the light-curve cannot be described as a sinusoid. In conclusion, it is evident that a 
significant variability is present in most of these data sets and that for Sycorax and Prospero the July 2005 data suggests 
the possibility to construct a periodical light curve. 

A period can be considered a good candidate if i.) x 2 (P) nas a local minimum; ii.) the amplitude A(P) of the associated 
sinusoid is significantly above the noise; iii.) the period P is not affected in a significant manner by aliasing with the 
sampling window. Of course one has to consider the fact that one period may be preferred to another one just by chance. 
Random fluctuations may lead to a different selection of the best fit period. We assess this problem by generating 
random realizations of time series with expectation given by the measured values of each sample and a fixed by their 
Gaussian errors. For each generation the period producing the minimum x 2 (P) has been determined. The probability 
of selection of each P, P se i(P), has been then derived. We then add a criterion iv) that a period P is selected as most 
likely if it has the maximum P se i(P). Tab. [3] reports the results of the fit. Note the difference in \ 2 between the best fit 
with a sinusoid and the x 2 m Tab. [U 



4 The Pfa (A > A^. Tit ) is the probability that random errors are responsible for the occurrence of a peak of squared amplitude 
A 2 > j4crit m the interval [P m i n , Pmax] of interest. 
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Before looking at the results, it has to be stressed that light-curves can be either single peaked or double peaked. In 
the second case the rotation period will be twice the light - curve period. We have not enough data to discriminate 
between these two cases, therefore the rotational periods of the observed objects could be twice the light-curve periods. 
Sycorax 

Fig. [5] on the left is the periodogram for Sycorax. The % 2 suggests that the three periods P « 3.6, 3.1 and 2.8 hr are 
favoured with a very high level of confidence (Pfa < 1CP 8 ). Bootstrap shows that the first period is the preferred one 
in about 96.6% of the simulations. The third period is chosen in less than 3.4% of the cases, the second one instead is 
chosen in less than 0.01% of the cases. For the best fit case we obtain also Bg — 21.676 ± 0.013, Vq — 20.849 ± 0.005, 
i?o = 20.276 ± 0.003 which are compatible with the weighted averages discussed in the next section and are fairly 
independent of P. The same holds for the "derotated colours" B - V = 0.828 ± 0.014, V - R = 0.573 ± 0.006. 
Prospero 

Fig. [8] shows on the right the periodogram for Prospero. Given a safe dependence between the phase angle and magnitude 
is not established for these bodies, we just used the data taken on July 2005 to evaluate the periodogram. Four periods 
are allowed at a 5cx confidence level (c.l.) for P ss 4.6, 3.8, 5.7 and and 3.3 hours, respectively. Bootstrap shows that the 
first peak is the preferred one in about 91% of the simulations. The second peak is chosen in less than 6% of the cases, 
the other peaks instead are chosen in less than 3% of the cases. Comparing the periodogram with the spectral window it 
is evident how the secondary peaks are close to alias of the diurnal 24 hr and 12 hr periods. Removing a 24 hr sinusoid 
from data before to perform the fit depresses the 24 hr peak, but not the 4.6 hr peak. On the contrary the removal of the 
4.6 hr component strongly depresses the spectrum in the range P k, 3 — 8 hr. By fitting separately the first and second 
night data, and avoiding the implicit 24 hr periodicity, the preferred period is 4.3 hr. Another way to filter the diurnal 
24 hr periodicity is to shift in time the lightcurve of the second night to overlap the lightcurve of the first. Even in this 
case periods between 3 and 5 hours look favoured. The spectral window for the data shows a leakage corresponding to 
P k, 4.3 hours. A secular variation may introduce power at periods longer than 48 hours which should leak power at 
P w 4.3 hours. The periodogram for simulated data with a linear time dependence has a peak in the P « 4.3 — 4.6 hours 
region, but in order to have an amplitude in the periodogram of 0.2 mag a peak-to-peak variation in the simulated data 
about seven or eight times larger than the peak-to-peak variation observed in real data is needed. Exclusion of B and/or 
I data, or fitting only the R data does not change significantly the results. The same holds if we remove the three R 
points with the largest errors. As a consequence, the data suggests P m 4.6 hr with A « 0.21 mag. The lower frame of 
Fig. [5] represents the variations of B, V, R and I magnitudes folded over the best fit sinusoid; For the best fit case we 
obtain also B = 24.584 ±0.123, V = 23.841 ±0.053, R = 23.202 ±0.020, h = 22.805 ±0.043 which are compatible with 
the weighted averages discussed in the next section and are fairly independent of P. The same holds for the "derotated 
colours" B -V = 0.743 ± 0.134, V - Rq = 0.639 ± 0.057, R Q - I = 0.397 ± 0.047. 
Setebos, Stephano and Trinculo 

For the other three bodies no strong evidences are found for a periodicity in the present data. However for Setebos the 
data may be considered suggestive of some periodicity with Pfa ~ 0.7%, the preferred period being P w 4.38 ±0.05 hours 
with A w 0.189 ±0.038 mag. For Stephano and Trinculo best fit periods are Pstcphano ~ 2 hours and P-ivinculo ~ 5.7 hours 
with amplitudes ^stephano ~ 0.459 mag and A-rrincuio ~ 0.422 hours but with Pfa of 22% and 44%, respectively. 

4. Averaged Magnitudes and Colours 

The derivation of colours would have to take in account the removal of rotational effects from magnitudes. Otherways 
systematic as large as peak-to-peak variations in the light-curve can be expected. Lacking a good light-curve we may 
apply two possible methods are: hierarchical determination of colours and compare weighted averages of magnitudes. 
The hierarchical method is based on the comparison of magnitudes from consecutive frames in the hypothesys that 
time differences are smaller than the light-curve period, so that rotational effects can be neglected. An example is the 
estimation of Setebos V-R by taking V3 and R15. In case in which one of the frames obtained with a given filter X is 
located between two frames of another filter Y, the Y magnitude at the epoch of which the X filter has been acquired 
can be derived by simple linear interpolation. An example is given by the estimate of V-R for Prospero in Nov 22, by 
interpolating R2 and R3 at the epoch of VI. Given in this way different estimates of each colours are obtained, weighted 
averages of such estimates are reported. The second method is based on the hypothesys that the light-curve has a 
periodical behaviour, and that it is so well sampled that weighting averages of magnitudes will cancell out the periodical 
time dependence. Of course the first method would be affected by larger random errors being based on a subset of the 
data. The second method is more prone to systematics. 

Colours derived by using the first method for Sycorax, Prospero and Setebos are presented in Tab. [5] When more 
independent estimates of the same colour are possible their weighted average is taken. 

Colours derived from the second method are listed in Tab. |4j The weighted averages of magnitudes for all of the 
satellites are listed in Tab. [1] In both cases, tables report just random errors and not the systematic calibration error, 
which amount to 0.018 — 0.022 mag, depending on the filter. 

It is evident how the results of the two methods are similar. Hence, we present for conciseness the results of the second 
one as commonly reported in literature. 
Sycorax 

We ob tain B — V = 0.89 ± 0.01 and V - R = 0.53 ± 0.01, which are both c ompatible within lrr with the iMaris et al~l 
(200ll) determinations but incompatible with iGrav. Holman fc Fraser I (|2004T ). 
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Prospero 

We have two sets of colour measures for Prospero, one from the July and the other from the November run. We obtain 
B — V = 0.89 ± 0.13, V - R = 0.63 ± 0.05 and R - I = 0.40 ± 0.04 from the July run, and B - V = 0.85 ± 0.13 , 
V-R = 0.59 ±0.07 and R-I = 0.33 ±0.10. The two set are consistent. With respect to lGrav. Holman k. Fraser I (|2004[ ). 
we obtain a redder V — R, and a compatible B — V. 
Setebos 

We obtain B — V = 0.74 ± 0.14 and V - R = 0.52 ± 0.04. In this case iGrav. Holman & Fraser I (|2004D obtained their 
data with the Keck II telescope and DEIMOS, which hosts a rather special filter set. Our B — V is compatible with 
their, while as for Prospero our V — R is redder. 

Stephano and Trinculo 

For these two extremely faint satellites, we could only derive the V — R colour, which is 0.73 ± 0.17 and 0.82 ± 0.43 for 
Steph ano and Trinculo, respectively. While our V — R for Trinculo is in good agreement with lGrav. Holman fc Fraser I 
(2004), their V — R for Stephano is very low, and inconsistent with our one. 

5. Summary and Conclusions 

In this paper we report accurate photometric B, V, R, I observations obtained with the VLT telescope in two consecutive 
nights in July 2005 of the Uranian irregular satellites Sycorax, Stephano, Trinculo and Setebos. Additional observations 
of Prospero obtained in November 22 and 25 of the same year are also reported. 

From the analysis of the data we conclude that S ycorax seems to displays a variability of m 0.07 mag, apparently 
larger than our previous result (|Maris et al. I . [200lT ). while the period of 3.6 hr is in agreement with our previous 2001 
determination, and the same is true for the colours, so it seems unlikely that the difference in amplitude can be ascribed 
to some systematic. If true, a possible explanation would be that in the two epochs two different parts of the same light 
curve have been sampled. But also it has to be noted that larger brightness variations have been not reporte d by other 
observers in the past years (jGladman et al. Lll998tlRettig. Walsh &; ConsoimagnoT 120011 : iRomon et al. 1 . 120011 ) ■ Prospero 
light-curve exhibits an apparent per iodicity of 4. 6 hr a nd an amplitude of 0.21 mag. The impact of such a sizeable 
amplitude is throughly discussed in (|Parisi et al. 1 . 120071 ) . Colours for Prospero obtained in July and November are in 
a quite good agreement, further assessing the goodness of the relative calibration. Setebos colours are only in marginal 
agreement with previous studies. In addition the Setebos light curve displays a significant variability but it is not possible 
to asses a good fit using a simple sinusoidal time dependence. Wether this is due to undersampling of the light-curve 
or to a non-sinusoidal time dependence can not be decided from these data alone. However, assuming a sinusoidal time 
dependence, our data are suggestive of a light-curve amplitude of rj 0.18 mag with a period of « 4.4 hours which will 
have to be confirmed or disproved by further observations. As for Stephano and Trinculo, the present data do not allow 
us to derive any sizeable time dependence, while the colours we derive are in marginal agreement with previous studies 
on the subject. 
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Fig. 1. Light curves in R for Sycorax. Data are for the nights of 2005, July 29 th (left) and 30 th (right). Squares in 
gray represents measurements of magnitudes for a common field star of similar magnitude. To avoid confusion error bars 
for the field star are not reported and the averaged magnitude is shifted. (See the electronic edition of the Journal for a 
colour version of this figure). 
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Fig. 2. Light curves R for Prospero. (See Fig. [JJfor comments). 
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Fig. 3. Light curves of Setebos. (See Fig. [JJ for comments). 
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Fig. 4. Light curves R for Stephano. (See Fig. [JJ for comments). 
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Fig. 5. Light curves R for Trinculo. (See Fig. [JJ for comments). 
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Fig. 6. Variability of the R light curves for irregular satellites (red circles) and field stars (light blue asterisks and 
gray band) as a function of their magnitude. Top frame for the peak-to-peak variation, bottom frame for RMS of the 
variation. (See the electronic edition of the Journal for a colour version of this figure). 
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Fig. 7. Time window for Sycorax (upper frame) and Prospero (lower frame). (See the electronic edition of the Journal 
for a colour version of this figure). 
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Fig. 8. LS test of periodicity for Sycorax (left) and Prospero (right) light-curve as a function of the period P. Upper 
frames display X 2 /Xmin> the re< ^ lme is the \ / Xmin f° r the case of a constant signal. The central frames are the estimated 
amplitude, a group of green lines denotes the amplitudes for Pfa — 10~ 4 , 10~ 5 , 10~ 6 , 10~ 8 and 10 -10 . In all the frames 
a black * marks the best fit point. Note that the diurnal 12 hrs and 24 hrs periods are excluded by the X 2 /Xmin despite 
the largest amplitude. The lower frames represents the overlap of the best reconstructed sinusoid with data ordered on 
a folded time scale, (t —r)/P, colours are blue for B, green for V, red for R and magenta for I. 
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Table 1. Log of observations. U.T. refers to the start times of exposures and are not corrected for light travel time. The Fit. 
column referes to the filter (BVRI) and frame number obtained with that filter (eg.: R3 is the third R frame for the given object 
in the serie). T cxp is the exposure time in seconds, the shortest exposures have been acquired to improve frame centering. 



Obj 


Epoch 


U.T. 


Fit. 


± cxp 

[sec] 


mag 


Sycorax 


Jul 29 


03 


27 


13 


Rl 


30 


20.282 ±0.015 


11 




04 


40 


37 


R2 


300 


20.270 ±0.005 


11 


// 


04 


46 


05 


R3 


300 


20.277 ±0.005 


11 




06 


46 


49 


VI 


300 


20.849 ±0.009 


11 




06 


52 


17 


V2 


300 


20.853 ±0.009 


11 




08 


29 


17 


R4 


300 


20.277 ±0.006 


11 


n 


08 


34 


45 


R5 


300 


20.313 ±0.007 


11 




09 


12 


19 


R6 


300 


20.325 ±0.008 


11 


n 


09 


17 


47 


R7 


300 


20.333 ±0.008 


Sycorax 


Jul 30 


02 


26 


25 


R8 


30 


20.214 ±0.023 


ii 


n 


06 


01 


03 


R9 


300 


20.283 ±0.006 


ii 


n 


06 


06 


29 


R10 


300 


20.289 ±0.006 


ii 


n 


06 


12 


15 


Bl 


300 


21.701 ±0.013 


ii 


ii 


08 


24 


01 


Rll 


300 


20.215 ±0.006 


ii 


n 


08 


29 


28 


R12 


300 


20.203 ±0.006 


ii 


n 


08 


35 


07 


V3 


300 


20.739 ±0.008 


ii 


n 


10 


09 


30 


R13 


300 


20.322 ±0.009 


ii 


n 


10 


14 


56 


R14 


300 


20.348 ±0.013 


Stephano 


Jul 29 


03 


34 


06 


Rl 


60 


23.502 ±0.183 


ii 


// 


03 


35 


34 


R2 


60 


23.960 ±0.297 


ii 


// 


04 


15 


10 


R3 


300 


25.419 ±0.262 


ii 


// 


04 


20 


37 


R4 


300 


24.414 ±0.242 


ii 


// 


04 


26 


04 


R5 


300 


24.271 ±0.210 


ii 


n 


04 


31 


31 


R6 


300 


24.177 ±0.195 


ii 


ii 


06 


16 


34 


VI 


600 


25.033 ±0.269 


Stephano 


Jul 29 


06 


27 


02 


V2 


600 


24.823 ±0.295 


ii 


// 


08 


03 


46 


R7 


300 


24.827 ±0.393 


ii 


n 


08 


09 


14 


R8 


300 


24.538 ±0.299 


ii 


ii 


08 


14 


41 


R9 


300 


24.211 ±0.207 


ii 


ii 


08 


20 


09 


R10 


300 


24.554 ±0.291 


Stephano 


Jul 30 


05 


35 


44 


Rll 


300 


24.190 ±0.166 


ii 


// 


05 


41 


11 


R12 


300 


24.119 ±0.157 


ii 


a 


05 


46 


39 


R13 


300 


24.056 ±0.148 


ii 


ii 


05 


52 


07 


R14 


300 


24.114±0.151 


ii 


ii 


07 


57 


58 


R15 


600 


24.302 ±0.196 


ii 


ii 


08 


08 


36 


V3 


600 


24.957 ±0.290 


ii 


ii 


09 


46 


41 


R16 


600 


24.285 ±0.269 


ii 


ii 


09 


57 


11 


R17 


600 


24.276 ±0.250 


Setebos 


Jul 29 


03 


05 


18 


Rl 


60 


23.230 ±0.176 


ii 


// 


03 


54 


07 


R2 


300 


23.274 ±0.088 


ii 


// 


03 


59 


35 


R3 


300 


23.337 ±0.092 


ii 


n 


04 


05 


03 


R4 


300 


23.383 ±0.092 


ii 


a 


05 


43 


04 


VI 


600 


23.796 ±0.080 


ii 


n 


05 


53 


32 


V2 


600 


23.768 ±0.083 


ii 


ii 


06 


03 


60 


V3 


600 


23.781 ±0.088 


ii 


ii 


07 


29 


43 


R5 


600 


23.196 ±0.080 


ii 


n 


07 


40 


12 


R6 


600 


23.242 ±0.088 


ii 


n 


07 


50 


40 


R7 


600 


23.336 ±0.099 


ii 


ii 


09 


25 


31 


R8 


600 


23.100 ±0.092 


ii 


ii 


09 


35 


58 


R9 


600 


22.847 ±0.074 




n 


09 


46 


25 


R10 


600 


22.903 ±0.077 


Setebos 


Jul 30 


02 


16 


30 


Rll 


60 


23.836 ±0.318 


ii 


// 


05 


01 


08 


R12 


600 


23.272 ±0.064 


ii 


n 


05 


11 


37 


R13 


600 


23.182 ±0.062 


Setebos 


Jul 30 


05 


22 


04 


R14 


600 


23.314 ±0.068 


ii 


// 


07 


02 


39 


R15 


600 


23.127 ±0.061 


ii 


a 


07 


13 


06 


R16 


600 


23.205 ±0.066 


ii 


a 


07 


23 


51 


Bl 


600 


24.455 ±0.132 
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Table 1. Continued. 



Obj 


Epoch 


U.T. 


Fit. 


[sec] 


mag 


11 


// 


07 


34 


30 


V4 


600 


23.649 ±0.081 


11 


// 


07 


44 


58 


V5 


600 


23.602 ±0.075 


11 


// 


09 


13 


23 


R17 


600 


23.223 ±0.077 


11 


// 


09 


23 


50 


R18 


600 


23.219 ±0.078 


11 


n 


09 


34 


18 


R19 


600 


23.249 ±0.084 


Trinculo 


Jul 29 


04 


15 


10 


R3 


300 


25.110 ±0.474 


ii 




04 


20 


37 


R4 


300 


24.862 ±0.378 


ii 


// 


04 


26 


04 


R5 


300 


24.544 ±0.272 


ii 


ii 


04 


31 


31 


R6 


300 


25.222 ±0.485 


ii 


ii 


06 


16 


34 


VI 


600 


25.925 ±0.671 


ii 


n 


06 


27 


02 


V2 


600 


26.189 ±0.901 


ii 


ii 


08 


03 


46 


R7 


300 


24.711 ±0.362 


ii 


ii 


08 


09 


14 


R8 


300 


25.237 ±0.568 


ii 


n 


08 


14 


41 


R9 


300 


24.626 ±0.330 


ii 


ii 


08 


20 


09 


R10 


300 


24.780 ±0.361 


Trinculo 


Jul 30 


05 


35 


44 


Rll 


300 


25.559 ±0.631 


ii 


a 


05 


41 


11 


R12 


300 


25.664 ±0.713 


ii 


ii 


05 


46 


39 


R13 


300 


24.996 ±0.386 


ii 


ii 


05 


52 


07 


R14 


300 


25.988 ±0.976 


ii 


ii 


07 


57 


58 


R15 


600 


24.550 ± 0.465 


ii 


n 


08 


08 


36 


V3 


600 


25.167 ±0.646 


Prospero 


Jul 29 


03 


29 


44 


Rl 


100 


23.275 ±0.130 


11 


ii 


03 


39 


16 


R2 


300 


23.175 ±0.082 


11 


ii 


03 


46 


44 


11 


300 


22.821 ±0.117 


11 


n 


05 


12 


00 


R3 


400 


23.005 ±0.057 


Prospero 


Jul 29 


05 


19 


09 


R4 


400 


23.125 ±0.064 


ii 


// 


05 


26 


27 


VI 


400 


23.680 ±0.079 


ii 


n 


05 


33 


34 


V2 


400 


23.763 ±0.085 


ii 


n 


06 


39 


14 


R5 


100 


23.239 ±0.120 


ii 


n 


06 


42 


09 


R6 


50 


23.160 ±0.150 


ii 


n 


06 


59 


28 


R7 


400 


23.383 ±0.102 


ii 


n 


07 


06 


36 


R8 


400 


23.362 ±0.100 


ii 


ii 


07 


13 


55 


12 


400 


22.996 ±0.140 


ii 


ii 


07 


21 


03 


13 


400 


22.966 ±0.144 


ii 


n 


08 


41 


56 


R9 


400 


23.087 ±0.084 


ii 


a 


08 


49 


04 


R10 


400 


22.979 ±0.072 


ii 


ii 


08 


56 


23 


14 


400 


22.552 ±0.090 


ii 


n 


09 


03 


30 


15 


400 


22.539 ±0.089 


Prospero 


Jul 30 


02 


13 


32 


Rll 


60 


23.779 ± 0.429 


ii 


a 


04 


18 


45 


R12 


98 


22.955 ±0.092 


ii 


n 


04 


21 


15 


R13 


210 


23.018 ±0.074 


ii 


n 


04 


25 


37 


R14 


20 


23.504 ±0.291 


ii 


ii 


04 


27 


11 


R15 


20 


22.929 ±0.182 


ii 


n 


04 


29 


12 


R16 


20 


22.949 ±0.205 


ii 


n 


04 


30 


18 


R17 


400 


23.205 ±0.075 


ii 


n 


04 


37 


26 


R18 


400 


23.296 ±0.080 


ii 


n 


04 


44 


43 


16 


400 


22.947 ±0.112 


ii 


ii 


04 


51 


51 


17 


400 


22.972 ±0.114 


ii 


ii 


06 


19 


32 


R19 


400 


23.411 ±0.080 


ii 


n 


06 


26 


39 


R20 


400 


23.392 ±0.079 


ii 


n 


06 


34 


20 


Bl 


300 


24.622 ±0.160 


ii 


n 


06 


39 


48 


B2 


300 


24.748 ±0.180 


ii 


ii 


06 


45 


40 


V3 


400 


23.867 ±0.092 


ii 


ii 


06 


52 


47 


V4 


400 


23.906 ±0.101 


Prospero 


Jul 30 


08 


42 


41 


R21 


400 


23.057 ±0.075 


ii 


// 


08 


49 


48 


R22 


400 


23.279 ±0.093 


ii 


n 


08 


57 


25 


18 


400 


22.792 ±0.113 


ii 


a 


09 


04 


32 


19 


400 


22.690 ±0.103 


Prospero 


Nov 22 


00 


44 


18 


Rl 


300 


23.193 ±0.077 


ii 




00 


50 


04 


Bl 


300 


24.583 ±0.181 
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Table 1. Continued. 



Obj 


Epoch 




U.T. 




Fit. 


[sec] 


mag 


11 


// 


00 


: 55 : 


49 


R2 


300 


23.240 ±0.081 


11 


n 


01 


: 01 : 


27 


VI 


300 


23.881 ±0.110 


11 


a 


01 


: 07 : 


06 


R3 


300 


23.310 ±0.090 


11 


n 


01 


: 12 : 


45 


11 


300 


22.726 ±0.157 


11 


a 


01 


: 15 : 


43 


12 


300 


22.925 ±0.243 


11 


a 


01 


: 18 : 


52 


R4 


300 


23.120 ±0.076 


11 


n 


01 


: 24 : 


38 


B2 


300 


24.710 ±0.218 


11 


a 


01 


: 30 : 


25 


R5 


300 


23.212 ±0.083 


Prospero 


Nov 25 


00 


: 36 : 


16 


V2 


300 


23.835 ±0.109 


11 


a 


00 


: 41 : 


54 


R6 


300 


23.194 ±0.080 


11 


a 


00 


: 47 : 


32 


13 


300 


22.838 ±0.184 


11 


n 


00 


: 50 : 


29 


14 


300 


22.969 ±0.207 


11 


a 


00 


: 53 : 


37 


R7 


300 


23.279 ±0.088 


11 


n 


00 


: 59 : 


23 


B3 


300 


24.686 ±0.203 


11 


n 


01 


: 05 : 


07 


R8 


300 


23.063 ±0.076 


11 


a 


01 


: 10 : 


45 


V3 


300 


23.837 ±0.108 


11 


n 


01 


: 16 : 


22 


R9 


300 


23.284 ±0.165 
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Table 2. Testing against random fluctuations. Columns 2 - 3 \ 2 an d significativity level (SL) assuming constant signal, 
4 - 5 x 2 an d SL assuming linear time dependence, 6-7 maximum amplitude and false allarm probability assuming 
periodical signal. Note that in all the cases a low probability (either SL or Pfa) denotes a high level of confidence in 
excluding noise fluctuations. 



Hypothesis HO 
Constant Linear Trend Periodical 



Object 


Xconst 


SL 


Xlincar 


SL 


^4 max 


Pfa 


Sycorax 


557 


< 1 x 10" y 


437 


< 1 x 10" a 


0.07 


< 1 x 10~ B 


Prospero 


92.7 


1 x 1(T 7 


85.2 


1 x KT 6 


0.22 


< 3 x 10~ 6 


Stephano 


45.3 


4 x 1(T 4 


44.9 


3 x KT 4 


0.36 


0.22 


Setebos 


59.2 


3 x 1(T 5 


59.2 


2 x KT 5 


0.19 


7 x 10~ 3 


Trinculo 


8.83 


8 x KT 1 


8.44 


8 x KT 1 


0.42 


0.44 



Table 3. Possible periods and amplitudes from fitting of Eq. [3] Column 1 x 2 f° r fitting, column 2 the confidence 
level (CL), Column 3 the best fit period in hours and its estimated internal uncertainty, Column 4 the corresponding 
amplitudes and their uncertainties. Solutions are ordered with increasing x 2 ■ 





Sycorax 




P 


A 


# 


x 2 


CL 


[hr] 


[mag] 


1 


90.734 




3.60 ± 0.02 


0.067 ± 0.004 


2 


134.88 


2a 


2.70 ±0.03 


0.065 ± 0.003 


3 


179.59 


3(7 


3.04 ± 0.02 


0.051 ±0.010 


4 


180.66 


3(7 


3.13 ±0.02 


0.051 ±0.010 




Prospero 




P 


A 


# 


x 2 


CL 


[hr] 


[mag] 


1 


27.288 




4.551 ± 0.040 


0.221 ±0.027 


2 


38.759 


2a 


3.827 ± 0.064 


0.201 ±0.029 


3 


53.286 


3a 


5.760 ±0.100 


0.162 ±0.090 


4 


55.079 


3a 


3.300 ±0.100 


0.121 ±0.090 



Table 4. Weighted averages of magnitudes 



Obj Run B 



Sycorax Jul. 

Prospero Jul. 

Prospero Nov. 

Setebos Jul. 

Stephano Jul. 

Trinculo Jul. 



21.701 ± 0.013 
24.678 ± 0.120 
24.651 ± 0.115 
24.455 ±0.132 



V 

20.807 ± 0.005 
23.788 ± 0.044 
23.851 ± 0.063 
23.713 ±0.036 
24.944 ± 0.164 
25.670 ± 0.413 



R 

20.276 ±0.002 
23.160 ±0.019 
23.196 ± 0.028 
23.192 ±0.018 
24.212 ±0.050 
24.855 ±0.117 



I 



22.760 ± 0.036 
22.864 ± 0.095 



Table 5. Colours derived from the hierarchical method 



Obj Run B-V V-R R-I 

Sycorax Jul. 0.915 ± 0.017 0.548 ± 0.006 
Prospero Jul. 0.793 ±0.138 0.590 ± 0.056 0.395 ± 0.046 
Prospero Nov. 0.792 ±0.171 0.635 ± 0.085 0.335 ±0.131 
Setebos Jul 0.831 ±0.143 0.444 ± 0.053 



Table 6. Colours derived from weigthed averages 



Obj 


Run 


B-V 


V-R 


R-I 


Sycorax 


Jul. 


0.893 ±0.014 


0.531 ±0.005 




Prospero 


Jul. 


0.890 ±0.127 


0.628 ± 0.048 


0.400 ± 0.041 


Prospero 


Nov. 


0.800 ±0.131 


0.655 ± 0.069 


0.332 ± 0.099 


Setebos 


Jul. 


0.741 ±0.137 


0.522 ± 0.041 




Stephano 


Jul. 




0.732 ±0.171 




Trinculo 


Jul. 




0.815 ±0.430 





